Entangling atoms in photonic crystals 
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We propose a method for entangling a system of two-level atoms in photonic crystals. The atoms 
are assumed to move in void regions of a photonic crystal. The interaction between the atoms is 
mediated either via a defect mode or via resonant dipole-dipole interaction. We show that these 
interactions can produce pure entangled atomic states. We analyze the problem with parameters 
typical for currently existing photonic crystals and Rydberg atoms. We show that the atoms can 
emerge from photonic crystals in entangled states. Depending on the linear dimensions of the crystal 
and on their velocity of the entangled atoms can be separated by tens of centimeters. 



I. INTRODUCTION 



Quantum entanglement is one of the most remarkable 
feature of quantum mechanics. Coherent control of the 
entanglement between quantum systems attracts lot of 
attention mainly because of its potential application in 
quantum information processing. Simultaneously exper- 
imental investigation of the entanglement allows us to 
test basic postulates of quantum mechanics and to an- 
swer fundamental epistemological questions. These ques- 
tions are related to the original Gedanken experiment of 
Einstein, Podolsky and Rosen jl| which triggered discus- 
sions about non-locality of quantum mechanics and mo- 
tivated experimental proposals to test whether quantum 
mechanics is the complete non-local theory. The first ex- 
perimental confirmation of the violation of Bell's inequal- 
ities H has been done with the help of entangled photons 
H . A weak point of experiments with photons is an insuf- 
ficient control of directions of emitted photons and small 
detectors efficiencies. This problem should be removed 
in proposals where highly excited (Rydberg) atoms are 
entangled. Probably the first proposal of such an exper- 
iment is described in Ref. 0. Other proposal have been 
presented in Re fs. ||. Authors of these schemes proposed 
techniques how to create entangled atoms in microwave 
single-mode cavities. Recently, controlled entanglement 
between atoms separated apr. by 10 mm interacting with 
an electromagnetic field in a high-Q cavity has been ex- 
perimentally realized . In addition, trapped ions have 
been created in entangled states fij] . 

In this paper we propose a simple scheme for entan- 
gling atoms in photonic crystals. We remind us that pho- 
tonic crystals are artificially created three-dimensional 
periodic dielectric materials which exhibit a frequency 
gap or several gaps in spectrum of propagating electro- 
magnetic (EM) waves |§|j]. An EM wave with its fre- 
quency from the gap can not propagate in the structure in 
any direction. Photonic crystals operating at microwave 
frequencies were successfully created in laboratories [ [To| . 
They consist of a solid dielectric and empty regions. The 
periodicity of a photonic crystal can be destroyed by re- 
moving or adding a piece of material which creates a 



defect EM mode in the structure. This mode is spatially 
localized around the region of the defect. The frequency 
of the mode and the spatial modulation of its electric field 
amplitude depends on properties of the defect Jll|-p^| . It 
means that one can adjust parameters of the defect mode 
by creating a suitable defect in the crystal. In particu- 
lar, the spatial dependence of the mode amplitude can 
be adjusted to particular needs. In quantum optics, de- 
fect modes in photonic crystals can be used similarly as 
high-Q single-mode cavities p5| , p!6[ . The quality factor 
of a single mode in a metallic cavity can be of order of 
10 8 or more and similar values can be reached for a sin- 
gle defect mode in a photonic crystal [^5) . Today three- 
dimensional photonic crystals are available only at mi- 
crowave frequencies. They can be used for experiments 
with Rydberg atoms, similarly as microwave cavities. 

In this paper we consider two interactions via which 
one can produce entangled atoms. Firstly, we show 
that at it is possible to generate entangled atoms with- 
out a defect mode, using the action of the resonant 
dipole-dipole interaction (RDDI) (ijjjl^] mediated by off- 
resonant modes of photonic band continua. Secondly, we 
explore the scheme in which the atoms become mutu- 
ally entangled due to the interaction with the defect-field 
mode. 

The paper is organized as follows: Basic features of the 
proposed setup are described in Section |n| In Section III 
we discuss how the atoms in photonic crystals can be 
entangled via the resonant dipole-dipole interaction. In 
Section [iy| we study in detail the entanglement of atoms 
which interact with a single defect mode in the photonic 
crystal. In Section [v| we conclude the paper with some 
remarks. 



II. SETUP OF THE SCHEME 

We consider two mechanisms via which a system of 
identical atoms can be entangled in photonic crystals. 
We assume that the atoms are modeled by two-level sys- 
tems having their transition frequencies in a photonic 
bandgap (PBG). 
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The first mechanism is the resonant dipole-dipole in- 
teraction (RDDI) mediated by off-resonant modes of the 
photonic-band continua [see the Hamiltonian (||)]. This 
interaction has been analyzed in detail by Kurizki [jisf] 
and John and Wang (l7]] as well as by John and Tran 
Quang fil}] . These authors have considered a system of 
two-level atoms. They have shown that if one of the 
atoms is excited and the other one is in its ground state, 
then they can exchange excitation in spite of the fact that 
their transitions frequencies are in a PBG and sponta- 
neous emission is nearly totally suppressed. The RDDI 
can be understood as an energy exchange via localized 
field |0. This light tunneling (or photon-hopping con- 
duction) can be very efficient when the distance between 
the atoms is much smaller than the light wavelength. 
The RDDI can occur either in a free space or in a cav- 
ity. However, in a free space the excitation is irreversibly 
radiated into the continuum of the field modes after a 
very short time (given by Fermi's Golden rule) and the 
entanglement between the atoms is deteriorated rapidly. 

The second mechanism is due to an excitation ex- 
change via a defect mode which is resonant (or nearly 
resonant) with the atoms. This type of interaction ex- 
plicitely involves a quantized defect mode and is de- 
scribed by the Hamiltonian (|l0|). 

These two interactions can also occur simultaneously. 
As we will see, the second mechanism is much more ef- 
ficient and allows a coherent control over the process of 
entanglement. The first mechanism can be neglected in 
many cases, especially when the atoms have their transi- 
tion frequencies near the center of a wide PBG and their 
distance is not much smaller than the wavelength of the 
resonant light. 

In what follows we describe the basic setup of the pro- 
posed experiment in the case when the atoms interact 
only via the defect mode. Let us assume that one of the 
three atoms (let say the atom A) is prepared initially in 
its excited state while the other two atoms (B and C) 
are initially in their ground states (see Fig.0). After the 
preparation the atoms are injected into cylindrical void 
regions of the crystal. We consider the photonic crystal of 
the geometry designed by Yablonovitch et al. JT^,|lJ| al- 
though other appropriate geometries can be used as well. 
The void cylinders intersect at the center of the crys- 
tal. The defect-field mode (located near the center of the 
crystal) is initially prepared in its vacuum state. Firstly 
the atoms propagate freely in the void cylinders outside 
the defect-field (this is due to the fact that the transition 
frequencies of the atoms lie inside the wide PBG). When 
the atoms enter the defect region they start to interact 
with the single defect-field mode. And then again, after 
they leave the defect region they evolve freely. If the ex- 
ited (ground) state of atom j (j — A, B, C) is denoted 
as \ej) (\gj)) and the n-photon state of the single mode 
defect field is denoted as \n) then the initial state of the 
system under consideration can written as 

|tf(0)> = \e A ) ® \g B ) ® \g c ) <S> |0) = \e A , g B , g c , 0). (1) 



When we assume that in the defect region the atom-field 
interaction is governed by the Hamiltonian in the dipole 
and rotating wave approximations (see below) then the 
final state of the system reads 

z 



It 




FIG. 1. A schematic description of the physical situation. 
We note that a particular geometry of the scheme and the 
type of the photonic crystal are not essential. The impor- 
tant feature of the crystal geometry is that it has straight 
"tunnels" so that atoms can traverse across the crystal. For 
concreteness, we choose the crystal with the geometry pro- 
posed by Yablonovitch et al. (references in Section Uf). We 
display only three of many cylindrical holes in the crystal. 
The cylindrical holes are drilled at the angle O = 35.26° with 
the vertical axis. We assume (except Section |il"]| ) that there 
is a defect of the crystal periodicity near the region where the 
holes are crossing. This defect is responsible for a single de- 
fect mode localized in the center of the crystal. The frequency 
of the defect mode lies inside a wide photonic bandgap of the 
crystal. Linear dimensions of the defect-mode region are com- 
parable with the lattice constant of the structure. The atoms 
A, B (and C, if needed) are injected into the three holes at the 
bottom side of the crystal at (approximately) the same time 
and with suitably adjusted velocities. The atoms are assumed 
to be two-level systems with their transition frequencies equal 
to the defect-mode frequency. We assume that transitions in- 
cluding other atomic levels can be neglected. The atom A 
is injected in its upper level |ca) and the atom B and C, in 
their lower levels |<?b) and \gc)- The initial state of the defect 
mode is vacuum state. The states of the atoms are detected 
at the exit from the crystal. We consider the following nu- 
merical values for the setup. The crystal is a cube of the 
side L w 20 cm. The frequencies of the defect mode and the 
atomic transition are Wo/(2tt) = w/(2tt) = 21.50651 GHz, i.e. 
the same as transitions used in experiments with microwave 
cavities. This frequency lies inside the wide photonic bandgap 
if the crystal is made from dielectric with refractive index (at 
microwaves) 3.6, volume filling fraction is 78% and the side 
of an elementary cube is a m 16.3 mm. We use this value of 
a to calculate the parameter k = ir/a [see Eq.(Ff2J)]. 
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= a{t)\e A ,9B,gc,0) + b(t)\g A ,e B ,gc,0) + c(t)\g A ,g B , e c , 0) + j{t)\g A , g B , gc, 1), 



(2) 



where i is the time at which we detect the internal states 
of the atoms at the exit of the crystal. The final values 
of the amplitudes a, 6, c and 7 depend on a particular 
setup of the experiment including the coupling param- 
eters and velocities of the atoms. For completeness of 
the description we specify trajectories rj(t) of the three 
atoms which can move along the axes of the three void 
regions 

r j (t)=r j (0)+v j t; j = A,B,C (3) 

with the vectors Tj(0) and Vj specified by their compo- 
nents as 

r A (0) = — {tan 6, -\/3 tan 6, -2} , 

v A = y j- sine, V3 sin 6,2 cos ej , (4) 

for the atom A. While for the other two atoms (B and 
C) we have 

r s (0) = j |tan9, V3tan9,-2|; 

w B = { — sin9, — V3sin 9, 2 cose| ; (5) 

and 

r C (0) = I {-tan 6,0,-1}; 

v c = vc {sin 6,0, cos 6} . (6) 

Here we assume the origin of the coordinates in the cen- 
ter of the cube crystal with the side of the length L; 6 
is the angle between the axes of the cylinders and the z 
direction. 



III. ENTANGLEMENT VIA RESONANT 
DIPOLE-DIPOLE INTERACTION 

In this Section we consider just two identical atoms 
(A and B) which move in the crystal as it is described 
above. Here we assume that there is no defect mode 
in the crystal. The atoms move inside the crystal with 
constant velocities. The recoil effect due to interaction 
with electromagnetic field is neglected because the atoms 
are relatively heavy particles. The interaction between 
the atoms and the electromagnetic field modes inside the 
crystal is described by the Hamiltonian in the electric- 
dipole approximation 

H = Tilu a 3 z + h u>xa\a\ 

3=A,B \ 

--^).DM-^(B).D(r B ), (7) 

where a\ and a x are the annihilation and creation oper- 
ators of the field mode labeled by A, D(r) is the trans- 
verse displacement-field operator and /J-{A) and n(B) are 
the atomic dipole operators. When the atomic transition 
frequencies are far from abrupt changes in the density of 
modes the Hamiltonian (Q) can be approximated as (for 
more details see Q) 

H cS = hLU ai + h(j AB a^a B + J BA a A af) , (8) 

J=A,B 

where a± are raising and lowering operators of the atoms 
(x = j4, B) and J AB is a matrix element for the effective 
description of the RDDI p7| . For qualitative estima- 
tions, we will use J AB evaluated under the assumption 
that density of electromagnetic modes is that of a free 
space. In this case we find (for more details see [p0[) 



fiJ AB = i4 e {A)ixf{B)- —[{5^ - 3RiRj)(cos k A R + k A Rsin k A R) - (5 V - R t Rj)k A R 2 cos(k A R)}, (9) 



where R is the distance between the atoms, k A = uj/c, 
fi e9 is the absolute value of the atomic dipole matrix el- 
ement, and Ri are the components of the unit vector 
starting at the position of the atom A and oriented to- 
wards the atom B. We assume summation over repeated 
indeces. We stress that the above expression for J AB is 
valid in a free space, but in the limit R <C A it can also 
be applied for photonic crystals [0, i.e. it can be also 
used for an order-of-magnitude description of the RDDI 
effects in photonic crystals. These effects are most impor- 
tant in the regime R <^ X when the free-space expression 
is valid also in photonic crystal. We will apply the Hamil- 



tonian (||) with J AB given by Eq.(^) also for description 
of propagation of the atoms in the crystal in the case 
when R > L. Even though the expression for J AB given 
by Eq.(^|) is not precise it provide us with a rather good 
picture of the RDDI effect. We note that in order to 
find more appropriate expression for J AB we would have 
to know the electromagnetic eigenmodes for the three- 
dimensionally periodic structure and the corresponding 
derivation of J AB is very complicated. 

In what follows we study the time evolution of the 
atoms initially prepared in the state l'J'(O)) = \e A ,g B ) 
which is governed by the effective Hamiltonian (|) with 
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time-dependent Jab (which is due to the fact that the 
atoms are moving through the crystal). We show that 
the RDDI can in principle be used for controlling the 
entanglement between atoms. We have solved the cor- 
responding Schrodinger equation numerically. We have 
used parameters typical for Rydberg atoms and currently 
existing photonic crystals. In Fig.|| we plot results for the 
time-dependent atomic populations. We have chosen the 
atomic trajectories similarly as it is specified in the pre- 
vious section but we added a small value (0.05 — 0.3 mm) 
to the initial xa(0) coordinate so that the trajectory of 
atom A is parallel but not identical with the axis of the 
cylinder. This prevents the collision of the atoms. The 
velocities of both atoms are 200 m s _1 . 



entanglement depends on the population of internal lev- 
els of the atoms and highest degree of entanglement is 
attained for \a(t)\ 2 = \b(t)\ 2 = 1/2. 

As seen from Fig. ||| the population of the excited state 
of the atom A depends on the minimal distance i? m ; n 
between the atoms during the passage through the crys- 
tal. From our numerical investigation it follows that the 
atoms are most entangled for i? m i n ~ 0.05 mm. However 
we note that with present techniques the controle over the 
position of atoms in the configuration considered here is 
about ±1 mm Q. Consequently, the RDDI is not very 
suitable for a coherent controle of entanglement between 
atoms in photonic crystals. In the following Section we 
consider entanglement via a defect mode when the cur- 
rently available precision control is sufficient. 
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FIG. 2. The time evolution of the population of the upper 
level of the atom A if the atoms interact according to effective 
Hamiltonian (^). The atoms move along the trajectories spec- 
ified in Section Three curves correspond to three different 
values of xa(0) specified by the minimal distance Rmin of the 
atoms during the passage. Both atomic dipoles are oriented 
in the x-direction. We evaluate the interaction only during 
the time interval when the atoms move near the center of the 
crystal in the cubic region of the side 2 cm. The atomic ve- 
locities are va — vb = 200 m s _1 and fi e3 je ~ 6.72 10 -7 m, 
where e is the proton charge. 

Taking into account that the physical conditions are 
chosen such that the electromagnetic field is adiabati- 
cally eliminated from the interaction [see the effective 
Hamiltonian (@)] the two atoms due to the unitarity 
of the evolution remain in a pure state \^(t))AB = 
a(t)\eA, 9b) + b(t)\gA, cb) with the amplitudes a(t) and 
b(t) which depend on the RDDI. From here it follows 
that due to the RDDI the two atoms become entan- 
gled. The degree of entanglement in the present case 
can be quantified with the help of the von Neumann en- 
tropy S = — Trfplnp] of each individual atom for which 
we have S = -\a(t)\ 2 In \a(t)\ 2 - \b(t)\ 2 In \b(t)\ 2 where 
| a(t) | 2 = 1 — |&(£)| 2 - In other words the degree of the 



IV. ENTANGLEMENT VIA A DEFECT MODE 

Let us consider the interaction of the atoms with a sin- 
gle defect-field mode in the dipole and the rotating-wave 
approximations. We assume that the distance between 
the atoms is always sufficiently large so that they do not 
interact via RDDI. The corresponding Hamiltonian can 
be written as 



H = huj <?i + Tiuj^a) a 

j=A,B,C 



ft [G(T j )a + a + G*(r j )aia 



(10) 



j=A,B,C 



where ujq is the mode frequency (which we assume to 
be equal to the atomic transition frequency u>), o± are 
atomic raising and lowering operators and and r b are 
the positions of the atoms. The position dependence of 
the coupling parameters G(tj) can be expressed as 



G(r j ) = G e-V j f( Tj 



fill 



where /(r) is the field-mode amplitude at the position 
r, e is the electric-field polarization direction of the de- 
fect mode and T>j is a unit vector in the direction of the 
atomic dipole matrix element of the atom j. It is known 
that the spatial dependence of a defect-mode amplitude 
is a function which oscillates and decays exponentially 
p"2[ . A particular profile of the spatial dependence of the 
defect mode can be adjusted via a properly generated 
defect of the periodicity. A rigorous calculation of the 
electromagnetic field in the presence of a defect in a 3D 
photonic crystal can be a difficult task. In this paper we 
use a model profile of the spatial dependence of the elec- 
tric field. Similar profiles have already been created in 
existing photonic crystals Q-Q. We note that for the 
purpose of the proposed experiment a complete informa- 
tion about the mode shape is not needed. The results 
of the experiment depend only on the shape along the 
trajectories of the atoms. In what follows we use the 
profile 
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/(r) = exp 



Rn 



-Rd. 



cf 



sin(k • r + $), 



(12) 



where Ro is the position around which the mode is local- 
ized, i?def is a parameter (defect-mode radius) describing 
the rate of the exponential decay of the mode envelope, 
<f> is a phase factor and k is the parameter describing 
spatial oscillations of the field mode. We chose its mag- 
nitude to be /c = n/a where a is the value of the side 
of an elementary cubic cell in the photonic crystal. Wc 
consider values of the constant i?def comparable with a. 
We estimate the value of the coupling constant Go from 
microcavity experiments ]16| 



G = 



n, 



(13) 



where V^av is the modal volume of the microcavity mode, 
V c s is the effective modal volume of the defect mode and 
SI is the vacuum Rabi frequency in the microwave exper- 



iment. The numerical values are [11 61: V cl 



11.5 cm 3 



and SI = 43 kHz. When we consider the transitions be- 
tween levels 63P 3 / 2 and 61D 3 / 2 of Rubidium atoms, then 
the atomic transition frequency is u/(2tt) = 21506.51 
MHz. Finally, the effective modal volume can be ap- 
proximated as 



V eS = -n(2R doi ) 3 . 



(14) 



Because the atoms are moving the coupling parameters 
depend on time [in what follows we will use the nota- 
tion Gj{t)] We consider positions of the atoms given by 
Eqs.([|) and (|J). In some cases we add a small value to 
xa (0) given by (Q) to prevent the atoms to collide in the 
center of the crystal. Details of the geometry of the pro- 
posed experiment are given in Section || and in Fig .|l|. 

Once we have specified all model parameters we can 
solve the Schrodinger equation for the system which is 
supposed to be initially prepared in the state [^(O)) = 
| eA , <7b j gc > 0) • Due to the fact that the number of excita- 
tions is an integral of motion in thepresent case the state 
vector at time t > has the form (||) and the correspond- 
ing Schrodinger equation can be rewritten into a set of a 
system of linear differential equations. These equations 
can be solved analytically for time-independent coupling 
constants Gj(t) which is not our case. Therefore we have 
to integrate the equations numerically. 



A. One atom 



We start our discussion with a problem when just a 
single atom (let say the atom A) passing through the 
crystal is considered. We assume that the atom is on 
resonance with the defect mode (i.e., u> = ujq). 
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FIG. 3. (a) The time dependence of the coupling Ga(£) between the defect mode and the single atom A when it moves 
along the axis of the cylinder with the velocity va = 500 ms -1 . The mode position and the geometry is given by the parameters 
Ro = 0, $ = rad and i?def = 10 mm. The parameter k = (0, 0, n/a) where a « 16.3 mm. The atomic dipole is oriented in 
the x direction (same as the field polarization). The integral ( |l7|) in this case is equal to zero. Consequently, the atom at the 
exit of the crystal is again in its initial state. We plot the time evolution of the population of the exited level of the atom in 
Fig.b. 



This corresponds to the Jaynes-Cummings model |2i| ] 
with a time-dependent coupling constant. The general 
solution of this model for real coupling parameter was 



found Sherman et al. p2fl . With the initial condition 
|4'(0)) = |e J 4, 0) the solution can be expressed as 
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= cos 



i sin 



G A {t')dt' 
G A {t')d£ 



\e A ,0) 



\9A,1) 



This implies for the atomic excitation 



G A (t')dt' 



(15) 



(16) 



In the case of the defect mode with linear dimensions 
much smaller than the side of the crystal we can use the 
approximation 



G A {t')dt' 



G A {t')dt'. 



(17) 



We note that this integral for a given choice of the profile 
function [see Eq.(|l2])] with the phase of the field mode 
$ = equals to zero. This means that the atom exits the 
crystal in the same state as it entered it. Obviously the 
defect mode also remains in its initial (vacuum) state. In 
Fig.||a we plot the time dependence of the coupling con- 
stant between the atom A and the defect mode. While 
in FigJ^b we present the time dependence of the corre- 
sponding excited-state probability. It is assumed that 
the defect is located at the center Ro = of the crystal. 
The other parameters are chosen such that $ = rad , 
k = (0,0, k), V A = e= (1,0,0) [see Eqs. ©,©]. We 
assume that the atom moves along the axis of the cylin- 
drical cavity. The velocity of the atom is chosen to be 
v A = 500ms _1 . From Fig. ||a we clearly see that the atom 
on its way through the crystal interacts with the defect 
mode just around the center of the crystal. The other 
important feature is seen from FigJ^b, i.e. The atom is 
transiently entangled with the defect mode in the center 
of the crystal. Nevertheless it leaves the crystal in a pure 
(unentangled) state. This effect of "spontaneous" dis- 
entanglement of the atom from the defect mode is very 
important when we consider creation of pure entangled 
state of two atoms. 



B. Two atoms 

Let us consider a situation when two atoms interact 
with the same defect mode as in the previous case. The 
atoms have their dipoles oriented along the direction e of 
the electric-field polarization. The velocity of the atom 
A is 500ms _1 . The time evolution of the corresponding 
atomic populations for various velocities of the atom B 
are plotted in Fig. ||. 

Firstly we consider both atoms to have the same veloc- 
ity (see Fig. ||a). In this case we assume that the atom A 
is displaced from the axis of the cylindrical hole through 
which it flies [i.e. we add 0.3 mm to x A (0) given by (^)] 
to avoid the influence of the RDDI between the atoms 
and their collision. We see that the atoms strongly inter- 
act with the field in the region of the defect. However, 



after the interaction the initial state of the system is ap- 
proximately restored (see the "stationary" values of the 
probability amplitudes a(r), b(r) and j(t) which are dis- 
played in the figures). It is interesting to compare Fig. [p3 
with Fig. [|a to see how the time evolution of the popu- 
lation of the atom A is modified by the presence of the 
additional atom B. We see that for the given set of pa- 
rameters the presence of the atom B does not influence 
the dynamics of the atom A significantly. 

Now we will study how the level population depends 
on the velocity of the atom B. From Fig. |] we see that 
for properly chosen velocity the interaction between the 
atoms mediated by the defect field can be pronounced. 
For instance, from Fig. ^b (here vb — 490ms~ 1 ) we see 
that not only the excitation of the atom B can be higher 
than the population of the atom A, but also the defect 
mode becomes partially excited and entangled with the 
atomic system. 

When the atom B has the velocity vb = 515ms -1 (see 
Fig. ||c) then the defect mode in the stationary limit is 
in the vacuum state [7(7") — — 0.0616i] and is (with high 
precision) completely disentangled from the atomic sys- 
tem. It is interesting to note that in this particular sit- 
uation the defect mode mediates transfer of most of the 
excitation from the atom A to the atom B. 

Let us assume now the velocity of the atom B to be 
Vb = 532.8ms -1 (see Fig. [|d). In this case the defect 
mode in the stationary limit is again in the vacuum state 
and is completely disentangled from the atomic system. 
Moreover the amplitudes a(r) and 6(r) are in this case al- 
most equal, which means that the atoms at the exit from 
the crystal are in the state I*) = (\e A ,gs) + \g A ,eB))/V2, 
i.e. they are prepared in a pure maximally entangled 
state. 

In the cases presented in Fig. || the phase factor $ of 
the defect mode is set to zero so that the integrals of the 
coupling constants G A (t) and Gs{t) over the trajectories 
of the atoms are equal to zero. The defect-mode radius 
Rdd — 10 mm. We have also studied the dynamics for 
other values of <E>, when the integrals of the coupling con- 
stants differ from zeros. In this case the disentanglement 
of the defect mode and the atoms is not so well pro- 
nounced, i.e. the defect mode becomes excited. We have 
also found a general feature: If the integrals of coupling 
constants are zeros and the coupling constants are small 
enough then the defect mode after the interaction is left 
in the vacuum state. However, if we increase the cou- 
plings (by decreasing the mode volume V e g) the defect 
mode can be left in an excited state [i.e. 7(7-) 7^ 0; see 
the expression for the state vector (||)]. Consequently, 
the atoms are left in a mixture state. 

We have also analyzed the situation when the defect 
mode is not located directly at the center of the crys- 
tal. Moreover we have assumed that <f> 7^ 0. It can be 
shown that even in this case it is possible to find a value 
we at which the atoms exit the crystal in a nearly pure 
maximally entangled state. 
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FIG. 4. The time evolution of the populations of excited levels of the atoms A (solid line) and B (dashed line). The 
atom A enters the crystal in the excite state, while the atom B is initially in the ground state. The atoms are injected into 
the crystal at the same time (see the setup presented by Fig.Q) with the velocity of the atom A being va = 500 ms _1 and 
Rdef = 10mm. The four plots corresponds to four different velocities of the atom B (their numerical values are shown in 
figures). In the case when the two atoms are assumed to have equal velocities (Fig a) we add a small value (0.3 mm) to 
the initial position xa(0) [see the expression (^)] to avoid the collision of the atoms in the center of the crystal. The other 
parameters are chosen same as in FigQ We write final values of the probability amplitudes into figures. From Fig. (d) we 
see that for properly chosen velocities of the atoms and times large enough (i.e. the atoms have already left the crystal) the 
defect mode is approximately in the vacuum state (it is disentangled from the atoms) and the atoms are prepared in a pure 
superposition state |*) ~ (|eA,5s) + \gA, e B ))/V^- 



C. Three atoms 

Let us consider the same setup as in our previous dis- 
cussion except we assume now three atoms flying through 
the crystal (see Fig. [I]). These three two-level Rydberg 
atoms (A, B and C) are injected into the holes at the 
bottom side of the crystal simultaneously. The atom A 
is initially in its upper level \ca) while atoms B and C 
are initially in their lower states \gs) and \gc)- The sin- 
gle defect mode is initially prepared in its vacuum state 
|0). The atoms move along the axes of the holes and 
interact with the defect mode in the central region of 
the crystal. The electric-field amplitude of the mode is 
given by Eq. (|T2|). We consider slightly asymmetric po- 
sition of the defect mode in the crystal (the reason is 



explained below). In Fig.|] we present plots of the final 
atomic populations versus velocities vb and vc while va 
is fixed at the value 500ms _1 . These plots show that 
adjusting atomic velocities we can obtain required prob- 
abilities such that in the final state (||) the probabil- 
ity amplitude 7(r) is equal to zero, which means that 
the defect mode is decoupled from the atomic system. 
The atoms are then in a pure superposition state. In 
particular, if we select velocities vb — 536.4 ms' 1 and 
vc = 527 A ms _1 , we obtain a final state with equal prob- 
abilities Ht)| 2 = |6(r)| 2 = \c(t)\ 2 m 0.33 (see Fig|). 
Square of |7(t)| gives the probability of the photon in 
the final state approximately equal to 0.02. It means 
that the atomic subsystem is in a good approximation 
decoupled from the field subsystem. 
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FIG. 5. Final atomic populations of the three atoms which interact with the defect mode in the crystal versus velocities vb, 
vc at fixed value va = 500ms _1 . The defect region is positioned unsymmetrically: Ro = (1,-3,2) mm. We chose the phase 
$ = and the other parameters (a and Rdei) same as in Fig.H. 



We have chosen an asymmetric position of the defect 
mode with respect to the center of the crystal because 
for the symmetric position we were able to obtain the 
"symmetric" result \a(r)\ 2 = \b(r)\ 2 = \c{t)\ 2 w 0.33 
only when two of the velocities are equal. In this case 
we face the problem of the collision of the atoms. We 
expect that a better choice of the defect geometry might 
produce a final state more disentangled from the field as 
is the case presented in Fig.|| 

We see from Figs.|^ that variations of the final atomic 
populations are rather robust with respect to changes in 
velocities, i.e. uncontrolled velocity fluctuations (which 
in experiments can be reduced up to 0.4 ms _1 do not 
deteriorate the predicted entanglement. 



V. CONCLUSIONS 

In this paper we have shown that atoms can be en- 
tangled in photonic crystals via dipole interaction medi- 
ated by off-resonant modes or via an interaction with a 
single defect mode. In the first mechanism (RDDI) the 



atoms can coherently exchange excitation while only a 
very small part of this energy is radiated into the field. 
However, this interaction might not be easy to control 
in an experiment because it requires a high precision po- 
sition control of the position of the atoms. The second 
mechanism (via a single resonant defect mode) is experi- 
mentally more promising because it can be realized with 
currently available microwave photonic crystals and with 
highly excited Rydberg atoms. 

We have shown that atoms can be prepared in pure en- 
tangled states and that the probability amplitudes of the 
generated superposition states of the atoms can be coher- 
ently controlled by varying the velocities of the atoms or 
by varying the orientations of the atomic dipole matrix 
elements. 

In our scheme of entanglement via defect modes in pho- 
tonic crystals the distance between the entangled atoms 
at the exit from media depends on the size of the media, 
the angle between the atomic trajectories, the atomic ve- 
locities and the life of the atoms. For the parameters 
used in this paper the distance between the entangled 
atoms is of the order of tens of centimeters. 
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Finally, we think that investigation of dynamics of Ry- 
dberg atoms in photonic crystals is an interesting com- 
plement to current experimental cavity quantum electro- 
dynamics. 
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FIG. 6. The time evolution of the atomic populations 
|a(i)| 2 (solid line), \b(t)\ 2 (long dashed line), and |c(t)| 2 
(short dashed line), when the three atoms are injected 
into the photonic crystal simultaneously with the velocities 
v A = 500 ms" 1 , v B = 536.4 ms^ 1 and v c = 527.4 ms~\ 
All other parameters are the same as in Fig|E|. We see that 
the three atoms exit the crystal in a state nearly disentangled 
from the field state - we plot the sum |o(t)| 2 + |&(i)| 2 + |c(t)| 2 
(dotted line) which is close to unity. 
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